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Abstract 

We calculate the squared gluon field strengths of a heavy q-q-pair in the model 
of the stochastic vacuum. We observe that with increasing separation a chromo- 
electric flux tube is built. The properties of the emerging flux tube are investigated. 
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1 Introduction 

As is well known, the behavior of QCD at large distances is most probably not con- 
trolled by perturbation theory, but genuine non-perturbative effects are supposed to 
play an essential role in that infrared domain. This is indicated as well by analytical 
considerations JT[ as by numerical calculations from lattice-QCD ||. It is of course 
the ultimate desire to obtain analytical results directly from the QCD-Lagrangian, but 
successes in that direction are limited. At the moment all analytical results are based 
on models. 

In the model of the stochastic vacuum (MSV) it is assumed that the complicated mea- 
sure which emerges from the QCD-action in the region of slowly fluctuating fields can 
be described by a well behaved stochastic process [|[ J| J|]. This process can be visu- 
alized as a fluctuating background field representing the " physical vacuum" . 
In the most general form of the MSV one assumes that the above mentioned back- 
ground field has a convergent cluster expansion, but for many specific applications 
phenomenologically useful results can only be obtained after further assumptions are 
made. One of these assumptions is that the stochastic vacuum can be described by a 
Gaussian process, i.e. only the 2-cumulant does not vanish. In a non-Abelian theory 
even the assumption of a Gaussian process has to be specified further, since the choice 
of the stochastic variables is ambiguous. 

Our paper is organized as follows: 

In section 2 we shortly present the model of the stochastic vacuum in a restricted form 
as it has been applied very successfully to soft hadron-hadron high-energy scattering 

i i- 

In section 3 we apply this restricted form of the model in order to calculate the squares 
of the field strengths of Euclidean Wegner- Wilson- loops. Especially we shall test the 
importance of some more technical assumptions described in section 2 and check the 
consistency of the approximations by investigating the Michael sum rule of the total 
energy stored in the field of the Wegner- Wilson- loop. 

In section 4 we shall present the results, discuss them and compare, where possible, 
with lattice gauge calculations. 

Throughout the paper we work in an Euclidean space-time continuum. 



2 The model of the stochastic vacuum 

A rather extensive description of the MSV as applied in this paper has been given in a 
recent publication [@|, so we quote here for convenience only the main points. 
The crucial ingredient of the MSV is the correlator between the gluon field strengths 
F^(x,w) parallel transported to a reference point w 

F^(x,w) = &~ l (x,w) F^(x) &(x,w), (1) 
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where ¥ fiU (x) is the Lie-algebra valued gluon field strength tensor and <&(x, w) the path 
ordered non-Abelian Schwinger string: 



V exp 



N*-l 
(7=1 



i<? / da (x — w) )1 A ll (w + <r(x — «;)) 



c 



iV c is the number of colors and unless otherwise stated we use the fundamental repre- 
sentation. 

We make the crucial approximation that the correlator < F^ u (x,w) F p B a (y,w) >a, 
i.e. the vacuum expectation value with respect to the non-perturbative background 
field, is independent of the reference point w. Then the most general form of the 
correlator is: 



§AB ^ q^FF > 

< g 2 F^(x,w) F*(y,w) > = ^ _ - y ^ \ {5^5 ua - 5^5 up )k D(z 2 /a 2 ) 



I d d \ 1 

+ - {-^-{zaSfip - ZpS^a) + —{z p 5 ua - z a 5 up ) \ (1 - k)Di{z 2 / a 2 ) I (3) 



where z = x — y. 

The fundamental parameters of the model are the correlation length a and the gluon 
condensate < g 2 FF > = < g 2 F p c u F^ v >. 

The factors in eq.(||) are chosen in such a way that D (0) = -Di(O) = 1. 

By evaluating the expectation value of a Wegner- Wilson-loop with the correlator (jU) 

and Gaussian factorization one obtains the potential for a static q-q-pair || 

nip n QQ N 

+ j^pdpj dT^-2KD(p 2 + t 2 ) + {1- ^D^p 2 + t 2 M. (4) 

In an Abelian gauge theory without monopoles, where the homogeneous Maxwell equa- 
tion must hold, only the second structure in eq.(H) can occur, i.e. we must have k = 0. 
However, in a non-Abelian theory there is no reason for k to be zero and it is the 
tensor structure proportional to K which leads to the linear rising potential (first term 
in equation and hence confinement. 

From comparison with hadron spectroscopy, high-energy scattering [@] and lattice gauge 
calculations || the input to the right hand side of eq.(||) is pretty well fixed: 

„. o . 9 , . f d 4 k k 2 ( k (z/a) 
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D ' (2 /°> - -^JwWTW^Vm) 
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A 4 



A 4 



(5) 



with A 4 = |-, A 4 = -8vr 2 , a = 0.35fm, < g 2 FF >= 2.39GeV 4 and k = 0.74. 



3tt 

The are the modified Bessel- functions. 




: D(z 2 /a 2 ) 
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Figure 1: The functions D(z 2 /a 2 ) and D\{z 2 jar) of eq.(5) 



3 The squared components of the gluon field tensor 



In this section we evaluate the tensor of the squared field strengths in the presence of 
a Wegner-Wilson-loop. 

We introduce as usual the Wegner-Wilson-loop over a contour C 



<W[C}>= trV < exp^-ig (j) A fl (x)dx fl ^ 



> 



(6) 



The plaquette P a p{x) is a Wegner-Wilson-loop over a square in a, [3— direction centered 
at x with side length Rp. 
One easily obtains 



P af i{x) = N c --R%g 2 Y,F^{x)F^{x) + 0(R%) 

c 



(7) 



where there is no summation over a and (3 on the r.h.s. 

The non-perturbative vacuum expectation value of P a p is related to the gluon conden- 
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sate: 

E < P «/3 > = 6N C - \r% <g 2 FF> + 0{R%) (8) 

Q</3 

To calculate the square of the gluon field strengths in the presence of a static quark- 
antiquark-pair we define the quantity [10] [)ll|] 



4 < W[C) P a/3 (x) > - < W[C ] >< P af3 (x) > 
R% < W[C] > 



Here C is the Wegner- Wilson-loop corresponding to a static q-q-pair with side length 

Rw in 3-direction and T in 4-direction centered at the origin (see figure . 

In the limit Rp — > the elements of the tensor (§) are the differences between the 

expectation values of the squared elements of the gluon field strength tensor in the 

presence of the q-q-pair and the fluctuations of the non-perturbative vacuum. 

In order to evaluate Ai 7 ^ (x) in the model of the stochastic vacuum we first transform 

as usual the path ordered line-integrals in the Wegner- Wilson-loops into surface ordered 



surface-integrals with the help of the non-Abelian Stokes-theorem [12 [. 



<W[C]> = tr V < exp l—ig <j> A fJi (x)dx^ > 

= trVs <exp(-ig J^F fiU (x,w)da^ > (10) 

The surface S is boardered by the loop C and contains the reference point w. The 
indices of the surface elements are restricted to /u < v throughout this paper. 
Expanding the exponentials in eq.(|9|) we obtain: 



surfaceordered 



^Fl (x) 



R P <W 



< n n F a l . . . J? a n n 2 pa pb > _ 

^ » jl\V\ ^ fl n V„ y [IV* pa ^ 



I ' I > ( 1 ) 



t a t b 



surfaceordered 



n=l 



E H)" / • • • / da% Vl ■ ■ ■ da™ Un tr [t« • • • f ] / / do£>£ X 
< ff"^ • • • F*l Vn > < g 2 F^F b pa > ) (11) 



2! 



t a t b 



Here we have neglected the higher terms in the expansion of P a p since they are in higher 
order of Rp. In eq.fllll) we have two loops and hence two surfaces each containing the 
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same point w. We have thus chosen the surfaces as sliding sides of a pyramid based 
either on the Wegner- Wilson-loop (da w ) or on the plaquette (da p ) (see figure^). 




Figure 2: The two pyramids touching each other at the reference point w 



We use the factorization of the color components of the field strength tensor as in [13| 



and 14], i.e. for example 



< F Cl (l)F C2 (2)F C3 (3)F C4 (4) > = < F Cl (l)F C2 (2) > < F c ' a {3)F Ca {A) > + 

< F Cl (l)F Ca (3) > < F C2 (2)F C4 (A) > + 

< F Cl {l)F Ci {4) > < F° 2 {2)F Ca {3) >, 
where F c >(i) = F^(x h w) (12) 

and analogously for all other expectation values. We obtain 

( surfaceordered 
EH) 2 " /•••/ <v<*Mt ai -t fls "]x 
n=l 



If 



da p „,do p { ~ iY 



pairs 



/w"-" pa 



(13) 



Of the 2n field strengths running over the pyramid of the loop C two are correlated with 
a .F^ on the pyramid of the plaquette and the other (2n — 2) among one another. As 
soon as the two field strengths of one of the latter (n — l)-pairs are not direct neighbors 
in surface ordering, this ordering is suppressed by at least a factor a/L as compared to 
the fully ordered ones (a is the correlation length and L the linear extension of the loop 
C). We now neglect in any order of the expansion of the exponentials all terms 
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which are not fully ordered, i.e. which are suppressed by at least a factor a/L. For all 
fully ordered expressions the traces in eq.(|l^) are identical and one obtains: 



Ai^3 (x) 



4 1 



E 



H) 2n + 2 • N t 



surfaceordcrcd 



R% < W > \^ 4(2AQ™(iV 2 - 1) 



Ml"! ' • • d(T ^nU2n X 



// 



ai,/3i---a„,/3 n 

E< a 2 F Cl F Cl > ■■■ (14) 

pairs 



< g 2 F^ n ~ 1 „ FV™" 1 „ fl > x < „ F^ >< g z F^ n +\ F °" +1 > 

3 ''"n-l" !.-! A'/3 n _l i/ /3„„i 3 f-an v a n V fj,0 n Vjj n pa 



TpC n TpCn 



2 pd+i pC„4 



where the iV c -f actors are due to the traces. 

Now the integrand is symmetric and we may replace the ordered surface integrations 
by unconstrained ones, correcting for the larger phase space by the factor l/(2n)!: 



Ai^3 (x) 



4 1 



E 



H) 2 "+ 2 • AT, 



#f> < W > \ ^ 4(2jV c )»(iV 2 - l)(2n)! 



• • • / dcr..,,., • • • da 



V < q 2 ^ 1 F Cl > • • • (15) 

pairs 



<r n 2 F Cn - 1 TpC n -i > v ^ n 2 F Cn F Cn ■> <r n 2 F c ' n + 1 F^+i •> 



Performing the integration over each correlator separately we obtain: 



Ai^ (x) 



E 



H) 2n + 2 • iV ( 



4 1 

^< W > \^ A{2N c ) n (N? - l)2«- 1 (n - 1)! 

2 



do^da™ < 9 2 F^ V F% > 



4H) 



1 



8-i^(iV 2 - 1) < W > 

frci 



2 



n-1 



da^dall<g 2 F^F%>) x 



\ ' ' ~^) | / / J^W r W ^ n 2 EiC pC . 

2 n n!(2iV c ) n ' ' ' ' ' /"-" /- T " 



2;; 



(16) 



In the approximation used the sum is just < ^ /> and the final result is 



Ai^ (x) 



#4 &N C {N£-1) 



< S^fS > 



(17) 
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We note that in the integral one surface element is connected to the Wegner- Wilson- 
loop (indicated by the index W) and one to the plaquette oriented in a, /3-direction 
(indicated by P). 




Figure 3: Parametrization of the sliding sides of a pyramid, s and r € [0, 1]. 
We now parametrize the surfaces as indicated in figure ||] and obtain with eq.(|3|): 



Ai^ (a) 



R% {N 2 - i)8N c 



^2 J dsw dsp dry/ drp 



surfaces 



< g 2 FF > 
12 



E 



p<cr 



d(s w ,T W ) 
(<W<W ~ S^Svp) nD(z 2 /a 2 ) + 



d(s P ,T P ) 



1 ( d d \ 1 

2 ( ~d7 <KZu ^ PbP ~ Zp5 ^ + ~dz~^ Zp ^ vc ~ z<7 ^ v p)) ~ K ) D i{z 2 /a 2 ) | 



(18) 

1 2 



Here z = £, w — £ p and the sum goes over the 16 different surface combinations 
parametrized by ^(sw, T w) and £u(s p ,t p ). 

Since we are interested in the limit R P — ► the integration over rp becomes trivial 
if we expand the expression around rp = 1/2. We found that in all cases the limit 
Rp — ► was practically reached for Rp < |o. 

Numerically the limit of static quarks was reached for T > 6a. But in order to study 
the static case we took the limit T — > oo analytically. 



4 Results and Discussion 



4.1 General results 



The reference point w plays a crucial role in the non-Abelian Stokes-theorem and it 
has to be common to both surfaces in order to apply the MSV. Its choice influences 
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the results. Even for the application of the MSV to calculate a simple Wegner- Wilson- 



loop the resulting area law refers to the area including the reference point |15[|. This 
is certainly one of the less pleasant features of the model and a consequence of the 
approximation of the complicated measure of the QCD-action by a Gaussian one. If 
we deform the surface, necessarily higher cumulants must be occur, since the result has 
to be independent of the surface. For the one loop case the MSV is in accordance with 
phenomenology and lattice calculations, if the minimal surface is chosen. We therefore 
have always chosen the minimal surface in the evaluation of a Wegner- Wilson-loop by 
the MSV in the Gaussian form. 

To calculate Ai 7 ^ (x) we proceed the same way. We choose the reference point in such 
a way, that the total resulting surface of the sliding sides of the pyramids is minimal. 
Since the surface of the plaquette is considered in the limit Rp — > 0, this means, that 
the reference point w has to be in the plane of the loop C. This will always be our 
standard choice. In order to investigate the influence of the position , we will later also 
study the case where the reference point has an arbitrary position on the connection 
line between the plaquette V and the loop. 

We first concentrate on the typical non-Abelian tensor structure determined by the 
correlation function D. 

It can be seen by symmetry arguments that for that case the difference of the square 
of the magnetic field strengths, i.e. the plaquette P a /3(x) with no time (x^) component, 
vanishes identically. This means that the magnetic background field is not affected by 
the static color charges. 

The electric field perpendicular to the Wegner- Wilson-loop is also practically not af- 
fected but only the difference of the squared electric field parallel to the loop. 
In both cases the squared electric field difference is negative, i.e. the presence of the 
static source diminishes the vacuum fluctuations. 

In figure 8 we display the value of the squared field strength parallel to the spatial 
loop extension AF| 4 (x) = — < g 2 E 2 {x^,r) > q q- V acuum as a function of the coordi- 



nates r = y xf + x 2 , and X3 (the results are rotational invariant around the X3-axes) for 
different spatial separations Rw (see figure [2] for the choice of coordinates) . 

The squared field strength reaches its saturation value AF34 (x) ~ 14j^f for a spatial 
extension of the Wegner- Wilson-loop Rw ~ 4a = 1.4fm. The transversal extension of 
the flux tube, defined by 



_ Jdrrr 2 AFf 4 (0,r) 
rMS = V JdrrAF 2 4 (0,r) ' (19) 



is practically independent of R\y and about 1.8 times the correlation length a (see 
figure |). 

As mentioned above, the quantity AF34 (x) is the difference between the squared field 
strength in the presence of a static q-q-pair and the non-perturbative vacuum field 
strength. In figure 10 we display the unsubtracted squared field tensor of the color 
source in units of the vacuum gluon condensate < g 2 FF >. 
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R w [a] 



Figure 9: The transversal extension tms of the squared electric field as a 
function of the quark separation Rw . 



As can be seen, the presence of the static source has only a moderate influence on the 
global vacuum-structure («10%). 

Now we discuss the second correlation function D\ . The contribution of this correlator 
to AF 2 g (x) is numerically much less important than that of the correlator D. For 
small distances this is mostly due to the factor (1 — k) = 0.26, but for large distances 
the fact that with this correlator no string is formed becomes decisive. 
Again the square of the magnetic field is not influenced by the static source, but now 
the squared parallel and perpendicular electric fields are of the same size and build up 
a spherical distribution. In figure 13 we show AF 2 4 (3:3, r) + AF| 4 (^3, r) + AF| 4 (x%, r) 
for the correlator D\. 



4.2 Michael sum rule 



The total field energy due to the presence of the static quark source is just the potential 
energy of the quarks. We thus have a relation between the expectation value of the 
Wegner- Wilson- loop and the spatial integral over the sum of Ai 7 ^ (x). This is one 
of the Michael sum rules [tig] . In our case we can not test this sum rule absolutely, 
since we do not know g 2 . We can however estimate the effective value of g 2 from 
potential models of quarkonia to g 2 ~ 6 [17]. The Michael sum rule yields than a very 
sensitive test for the some more technical approximations made to calculate Ai 7 ^ ( 



x). 



We remind especially of the special choice of factorization and the neglect of all terms 
suppressed by at least a factor a/L in the expansion of the Wegner- Wilson-loop. 
In figure |l4| we display the numerically integrated energy density ■^{/S.F 2 i + A.F24 + 
AF 3 2 4 ) with g 2 = 7.2 and compare it to the q-q-potential (H) as determined from the 
model of the stochastic vacuum. 



We observe an almost perfect agreement between the two determinations even for 
distances which are small as compared to the correlation length a. (The error of the 
points is about the same size as the points due to the numerical integration.) The 
agreement in form and the very reasonable value of a s ~ 0.57 for the effective strong 
interaction coupling show, that the approximations of the MSV are well consistent. 
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Figure 14: The total energy stored in the field (dots) compared with the 
potential of a q-q-pair as obtained by eq.(4) 

4.3 Color-Coulomb contribution 

Having determined the effective a s we are in the position to obtain an approximation to 
the full QCD-expression by adding to the non-perturbative results the color-Coulomb 
contribution. The result for the full energy density is displayed in figure 17 for Rw = 2a 
and Rw = 9a. 



For R\y = 0.7fm the perturbative contribution increases the energy density in the 
middle of the flux tube from 0.65^^- to 4.2^^ and even for Rw = 1.4fm this increase 
goes from I-Ite^t to 1-9-g^r. Since lattice calculations are typically working with 
Wegner- Wilson-loops of spatial expansion Rw ~ lfm the measured field strength is to 
a large extent a perturbative one. The magnitude of our full energy densities (SU(3)) 
agrees qualitatively with the scaled lattice results for SU(2) |l]| [19] ]po|l . 



4.4 Dependence of the results on the reference point 



In a non-Abelian gauge theory it is essential to transport the color content of the field 
strength correlator to a fixed reference point in order to obtain gauge covariant results. 
Such a reference point also occurs quite naturally in the non-Abelian Stokes-theorem 
as starting point of the surface ordering [15]. The final result for the expectation value 
of a Wegner- Wilson-loop must be independent of the choice of the reference point, but 
the result of the MSV calculation does depend on it, if the reference point is moved 
outside the surface under consideration. As discussed before the minimal surface is 
most adequate for the simplification of the measure of the QCD-action to a Gaussian 
one. Thus this has been our standard choice. Now we want to test the dependence of 
the results on the reference point, varying its position on the connection line between 
the plaquette and the loop as indicated in figure 0. a=0 corresponds to our standard 
choice of w\ on the loop and for a=l t he reference point is inside the plaquette. 



In figure |l9] we display the squared electric field parallel to the quark axis in the 
transversal plane between the two quarks (A-F34 (0, r)). In figure ^ we show the same 
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Figure 18: Definition of the parameter a which describes the position of the 
reference point w. a=0 corresponds to our standard choice and for q=1 the 
reference point is inside the plaquette. 



field as a function of the position between the quarks with a distance of one correlation 
length from the q-q-axis (AF| 4 (#3, la)). For the parameter a we took the values 
a = i i 2 1 




Figure 19: The squared electric field parallel to the q-q-axis in the transversal 
plane between the two quarks (AF34 (0, r)) for a quark separation Rw = 6a 
and for different positions of the reference point w (see figure 18). The solid 
line is for a = and the values a = \ , \ , | , 1 are plotted by subsequent shorter 
dashes. 
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Figure 20: The squared electric field parallel to the q-q-axis as a function of 
the position between the quarks with a distance of one correlation length from 
the q-q-axis (AF34 (x%, la)) for a quark separation Ry/ = 6a and for different 
positions of the reference point w (see figure 18). 

As can be seen, the results are rather stable for a < | and the variation of a has only 
a limited numerical effect. 



4.5 Summary 

The MSV in its simplest form makes definite predictions on the non-perturbative con- 
tributions of the chromodynamic fields of a static q-q-pair. The first result is that 
the expectation values of the squared magnetic and perpendicular electric field are not 
affected by the quark pair whereas the squared electric field parallel to the quark axis 
is diminished. This is in qualitative agreement with lattice gauge calculations of the 



connected linear field strengths pif. The vacuum structure is quite stable, because the 
fluctuations of < g 2 FF > are decreased by only 10%. 

With increasing quark separation a flux tube with constant width (tms ~ 0.64fm) and 
energy density of about lj^j- is built. 

The approximations in the MSV made here respect the Michael sum rule with an 
effective strong coupling as=0.57. By adding the perturbative contribution with this 
coupling to the non-perturbative one we obtained the total energy density in agreement 
with lattice SU(2) gauge calculations U] II M\- 
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Figure 4: Rw = la, a =0.35fm Figure 5: Rw = 4a, a =0.35fm 
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Figure 10: The unsubtracted squared field tensor < g~FF > q q in units of 
< g 2 FF >vacuum for a quark separation Rw = 9a. 



Figure 11: Ryy = 7a, a =0.35fm Figure 12: R\y = 9a, a =0.35fm 

Figure 13: Difference of the squared electric field (— < g^E^(x^,r) >qq-vac) 
in caused by the correlator D\ for two different quark separations Rw- 




Figure 15: Rw = 2a, a =0.35fm Figure 16: Rw = 9a, a =0.35fm 

Figure 17: The full energy density of a q-q-pair caused by the non-perturbative 
correlator and the color-Coulomb contribution with as = 0.57 in for two 
different quark separations. 



